Coreless vortex ground state of the rotating spinor condensate 
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We study the ground state of the rotating spinor condensate and show that for slow rotation the 
ground state of the ferromagnetic spinor condensate is a coreless vortex. While coreless vortex is 
not topologically stable, we show that there is an energetic threshold for the creation of a coreless 
vortex. This threshold corresponds to a critical rotation frequency that vanishes as the system size 
increases. Also, we demonstrate the dramatically different behavior of the spinor condensate with 
anti-ferromagnetic interactions. For anti-ferromagnetic spinor condensate the angular momentum 
as a function of rotation frequency exhibits the familiar discrete staircase behavior, but in contrast 
to an ordinary condensate the first step is to the state with angular momentum 1/2 per particle. 
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I. INTRODUCTION 

The quantization of circulation for a Bosc-Einstcin 
condensate (BEC) is well known jij. This quantization 
implies a crucial role for quantized vortices as carriers of 
angular momenta. If the condensate rotates with angular 
frequency £1, the ground state of the system is a config- 
uration of quantized vortices and the number of vortices 
increases with the rotation frequency Energeti- 
cally a state with a single vortex becomes favored only 
after exceeds certain threshold £l c that depends on the 
system parameters in a fairly complicated manner j7j . 

While these properties are well studied, the creation 
of spinor condensates |||| has added an important ex- 
tra twist. In a spinor (or vectorial) condensate with spin 
S the spin-degree of freedom is not frozen and the con- 
densate order parameter is a vector £ with 2S + 1 com- 
ponents. This freedom with the spin-degree has some 
intriguing consequences. For example, in inhomogeneous 
magnetic field this vectorial nature of the order param- 
eter can even imply vortex ground states JTc( ] . Crucially 
for our purposes, such circulation of an order parameter 
does not have to be quantized jll|] . If circulation is not 
quantized angular momenta of the system does not have 
to be carried by quantized vortices. 

In this paper we address the question of the ground 
state of the rotating spinor condensate. We show that 
the ground state of the ferromagnetic spinor condensate 
at small rotation frequencies is a coreless (non-singular) 
vortex and we analyze its properties. Coreless vortices 
have also been studied by Mizushima et al. jljj who in- 
vestigated rotating spinor condensates with a fixed mag- 
netization. We mainly concentrate on the pure rotational 
ground state with no external magnetic field. Further- 
more, we study the rotating spinor condensate with po- 
lar (anti-ferromagnetic) interactions and show how sig- 
nificantly the behavior of the system changes when only 
the sign of the antisymmetric interaction parameter is 
altered. In both cases there is a critical rotation fre- 



quency above which the ground state carries non-zero 
angular momentum. However, in the ferromagnetic case 
the amount of angular momentum can be infinitesimal 
whereas the corresponding polar condensate has a quan- 
tized angular momentum (in trap units) 1/2 per particle. 
The vortex structure of a polar spinor BEC has been pre- 
viously studied by Yip jl3| focusing for a single constant 
rotation frequency. 



II. MODEL 

The Hamiltonian for the three-dimensional spin-1 con- 
densate in a frame rotating with frequency f2 is [|l2 14 1(| 



H= I d3r w « ' Wi + y W 
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where m is the atomic mass, u)t r and lo z are the trap- 
ping frequencies, L is the angular momentum operator, 
and F Q is the angular momentum matrix. The param- 
eters A s and A a are system specific parameters that de- 
pend on the two different scattering lengths and a : 

x = W(a +2a 2 ) ^ = jgft^sa Tn mean 

field limit our aim is to find the fields ^,(r) = ipi(r) 
that minimize this Hamiltonian. Minimization leads to 
the three-component Gross-Pitaevskii equations jl5) . We 
choose lu z = (2tt)1480 Hz and uj tr = (2tt)62 Hz. When 
the trap is filled with 1.7 x 10 4 Rb 87 atoms the condi- 
tion n < hto z is satisfied. Consequently, in z-direction 
we assume the system to be in the ground state of a 
parabolic trap with no dynamics and our numerical cal- 
culation becomes effectively two-dimensional. This as- 
sumption scales the scattering lengths by the constant 



1 



n 



For convenience we choose the unit of 



length as L = 



and the unit of time as tq = 1 / iot r 



It is well known that the ground state structure ib_= 
(ipi , ipo , 1 ) T = \fn ( depends on the sign of the A a |0 . 
In the absence of a magnetic field, if A a > the energy 
is minimized (with a suitable choice of axis) with a non- 
magnetized spinor £ = (0,1, 0) T . This state is referred 
to as the polar state. If A a < 0, a fully magnetized fer- 
romagnetic state C = (1) 0, 0) T is favored. In the absence 
of external magnetic field we can freely rotate the order 
parameter (globally) without any physical changes. In 
particular we can rotate the ferromagnetic ground state 
and get an identical order parameter expressed in terms 
of the Euler angles (a, 0, r) of the rotation 
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The superfluid velocity is defined by v s 
we make the rotation local and choose r 
4> = tan~ 1 (x/y), the velocity field of the condensate is 
given by v s = ±(1 - cos0)$ (l|. If 0(r) = the ve- 
locity field vanishes, but as 0{r) increases the velocity 
field starts to resemble that of an ordinary singular vor- 
tex. Such a structure is called a coreless vortex which 
was first predicted by two studies for superfluid 3 He. For 
the Mermin-Ho (MH) |ll| vortex the bending angle 
must be ? at the boundary of the condensate and for 
the Anderson- Toulouse (AT) jlT]] vortex must be ir. 
For the MH vortex in liquid 3 He the boundary value is 
imposed by the walls of the container |Q . 



III. FERROMAGNETIC GROUND STATE 

The angular momentum per particle in a condensate 
with a coreless vortex can be anything between and 2. 
Also, a coreless vortex can be continuously transformed 
into a state with vanishing angular momentum. This 
happens by simply transforming the function 0{r) into 
zero. Due to these reasons it seems plausible that when 
A a < 0, a coreless vortex might play a role as a carrier of 
angular momenta at small values of O. We have verified 
this reasonable conjecture by solving the relevant two- 
dimensional GP equations numerically for 1.7 x 10 4 Rb 87 
atoms without any additional assumptions: 



iti 



iti 
if i 



dt 
dip 



= + A a + 1^-1 + |Vo|V 



dt 
dt 



where 

C = + V trap ~nL+X s (l^-ll 2 + |V^0| 2 + IV'll 2 )- 

The lowest energy configuration at small values of SI is 
a coreless vortex (see an example of the superfluid ve- 
locity field in Fig. [j]). A priori, there is no reason to 
expect a critical rotation frequency for the existence of 
a coreless vortex. Nevertheless, by solving the GP equa- 
tions we saw that a critical frequency Q c ~ 0.06 had to 
be exceeded before a finite sized coreless vortex became 
energetically favored. Below this critical frequency en- 
ergy was minimized with the function 0(r) being zero 
everywhere. 
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FIG. 1. Superfluid velocity field inside a rotating ferromag- 
netic spinor condensate for the rotation frequency Q = 0.09. 
The angular momentum per particle is L/N ~ 0.6. 

Assuming the existence of a coreless vortex the Hamil- 
tonian in trap units is 



H C v [n, 0) = / d 2 rn 



1 2 
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where = Q/u> tr . When the rotation frequency is low 
and the density is large, it is reasonable to assume a 
density distribution almost independent of the function 
and that this density distribution has the familiar in- 
verted parabola shape of the Thomas-Fermi approxima- 
tion fl9[| . The 0{r) minimizing the energy functional (^) 
must satisfy the following differential equation 



£Vo + A a (2^_iV5 + |V-i|Vo + l^i |Vo) (3) 
A/>1 + Aa (^oCl + hM'tyi + IV'ol 2 ^ 
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This equation is solved numerically to get the lowest en- 
ergy configuration. Some results are shown in Figure |^. 
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By using the same parameters as in the numerical sim- 
ulations the energetical studies give ft c — 0.058 for the 
critical rotating frequency. This result is extremely close 
to the corresponding value from the GP equations. Also, 
the values of angular momenta as a function of f2 are 
nearly identical as can be seen in Figure |^. 

To gain more insight, Eq. can also be studied vari- 
ationally. With a trial (3(r) — r/L the resulting integrals 
can be solved. While the resulting expression is quite 
complicated, it is remarkable that such a simple varia- 
tional model gives results that are qualitatively the same 
as the results based on the Eq. (||). It turns out that the 
the existence of an energy minimum depends (in trap 
units) on the magnitude of the number A = R^ F fl. En- 
ergy mimimum only exists if A > 3.05 ± 0.05. Thus the 
critical frequency predicted by the variational model with 
the parameters used before is Q c = 0.064±0.001, i.e. very 
close to the value calculated with more accurate meth- 
ods. As the critical frequency is inversely proportional 
to the square of the Thomas-Fermi radius we see that 
the non-zero critical rotation frequency is a finite size ef- 
fect. This property is also confirmed by the numerical 
solutions of the GP equations. 




FIG. 2. The bending angle (5{r) for three different rotating 
frequencies: Q. = 0.07 (solid line), Cl — 0.08 (dashed line) and 
Q — 0.10 (dashdotted line). The corresponding values for 
angular momenta are L/N = 0.24, L/N = 0.42, L/N = 0.74. 
In this case the Thomas-Fermi radius is r = 6.9. In order to 
prevent the divergence of at this point the second term in 
the right hand side of Eq. (0) is put to zero at r — 0.99Rtf- 



A. Ground state as a function of fi 

As the rotation frequency is increased, the slope of 
(3 increases initially (to a good accuracy) linearly with 
the rotation frequency. This can be easily understood 
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FIG. 3. The angular momentum of a ferromagnetic spin-1 
condensate obtained by solving the GP equation (dotted line 
with plus signs) and by solving the differential equation (^) 
(solid line). 



since initially a larger slope implies larger populations 
for the state with the highest winding number, namely 
m = — 1 state. But if the slope becomes too large 0(r) 
will be larger than tt inside the condensate. If this hap- 
pens, populations at states with lower winding numbers, 
namely m = 1 and m = states, start to increase. This 
implies that angular momentum will decrease. Physically 
this is nonsense since it is clear that angular momentum 
per particle must increase monotonically with the rota- 
tion frequency. Consequently we expect an upper critical 
frequency O up when the coreless vortex ceases to be the 
ground state of the system. 

We can get a rough estimate for this upper limit by 
assuming that the slope depends linearly on f2 and ex- 
trapolating the behavior to higher rotation frequencies. 
For the parameters used before this gives us an upper 
critical frequency Q up ~ 0.15, which is again in good 
agreement with the solutions of the GP equations j2(J . 

The ground state at higher rotation frequency is more 
complicated. The superfluid velocity field slightly above 
the upper critical rotation frequency is still rotationally 
symmetric, but the order parameter is no longer a simple 
mixture of components with winding numbers 0, 1, and 
2. In Figure ^ we show an example of the order parame- 
ter when Q — 0.17. In the m = component we can see 
a regular array of 4 vortices. The m = ±1 components 
fill the vortex cores. Although it is not clear from the 
figure, the m = ±1 components have a nonzero angular 
momenta. The velocity field of this configuration is ro- 
tationally symmetric and qualitatively resembles that in 
Figure [l] 

In reality f3(r) is not linear, but its slope vanishes close 
to the edge of the condensate. It seems that the maxi- 
mum boundary value for the bending angle is (3 max = | ir 
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FIG. 4. Ground state structure of the ferromagnetic spinor 
condensate when Ct — 0.17. The parameters used to create 
these figures are the same as for Figure In the top most 
figure we show the density of the m = 1 component, in the 
middle the density of the m — component, and in the lowest 
the density of the m = — 1 component. Light color indicates 
high density. Total angular momentum per particle of this 
configuration is about 1.85. 

and AT-vortex is never the ground state of the system. 
We believe that this, somewhat surprising, boundary 
value is caused by a subtle effect due to kinetic energy. 
Ordinarily kinetic energy of the BEC can be safely ig- 
nored, but close to the boundary of the condensate where 
gradients are large its effects are pronounced. We conjec- 
ture that the coreless vortex structure with j3 ma x > §tt is 
dynamically unstable and that this instability is due to 
the physics close to the condensate edge. But we freely 
admit that at the time of writing this paper, a detailed 
understanding of the spinor condensate dynamics close 
to the boundary of the condensate is still lacking. 

IV. ANTIFERROMAGNETIC GROUND STATE 

It should be understood that the results change dra- 
matically with the sign of A a fL3|. While a coreless 
vortex can carry an infinitesimal amount of angular mo- 
mentum, this state is not energetically favored if A a > 0. 
For polar condensates the ordinary non-rotating ground 
state C = (0, 1, 0) T remains the ground state until a cer- 
tain critical frequency is reached. This critical frequency 
is somewhat higher than the corresponding critical fre- 
quency for a (ferromagnetic) coreless vortex, but well 
below the critical frequency for the stabilization of an 
ordinary singly quantized vortex. 

Just above this threshold the ground state of the sys- 
tem is a state with angular momentum 1/2 per particle 
(see Figure ||) . This state is locally magnetized and can- 
not be represented as a local rotation of the non-rotating 
ground state. This is easily understood, since there is 
no such local rotation of the polar state that results 
in a texture having non-zero angular momenta. While 
the state is locally magnetized, the number of atoms at 
rn = ±1 states is still the same and the system is glob- 
ally non-magnetized. We can choose a convenient basis 
by a global rotation of the numerically calculated wave- 
function. Such a procedure reveals that this state can 
be intepreted as a state where population is evenly split 
between m = ±1 states. Another one of these states 
carries a vortex and the other one fills the vortex core. 
Incidentally, the dynamical stability of this structure was 
predicted in Ref. |21j . Here we have shown that within a 
certain range of rotation frequencies, this state is also the 
minimum energy configuration. Very interesting point is 
that the winding number 1/2 vortex of the polar spinor 
condensate can viewed as the BEC counterpart to the 
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FIG. 5. The superfluid velocity field inside a rotating 
spinor condensate for two different rotation frequencies. Scat- 
tering lengths are such that the non-magnetic (polar) struc- 
ture is favored and Thomas-Fermi radius is Rtf ~ 6.9. In the 
upper figure the rotation frequency is Q. — 0.09 and angular 
momentum per particle is 1/2. In the lower figure Cl — 0.1 
and angular momentum is close to 0.8. 

Alice string familiar from the particle physics p2| . 

Even though the behavior of the system depends 
strongly on the sign of A a , it does not depend strongly on 
its magnitude. A a could be increased threefold without 
large qualitative or quantitative changes. Also, the num- 
ber of particles could be changed by, at least, a factor of 
two without important qualitative changes. The range of 
rotation frequencies when this state is the ground state 
is quite narrow. In our simulations this range was only 
on the order of 0.01. If the rotation frequency is beoynd 
this narrow range, the ground state is a state where pop- 
ulation is evenly split between m — ±1 states, but in 
this new state both components carry a vortex. These 
vortices are not overlapping and the total angular mo- 
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mentum per particle is less than 1. These results are 
analogous to the results for scalar condensates, where it 
has been shown that in a parabolic trap a multiply quan- 
tized vortex is never the ground state |||| . In a parabolic 
trap it is always energetically favorable to break a mul- 
tiply quantized vortex into pairs of singly quantized vor- 
tices. 



V. MAGNETIC FIELD 

So far we have assumed that the chemical potentials of 
the different m-states are equal. In this section we discuss 
the effect of an external magnetic field pointing along the 
rotation axis. If the magnetic field has components also 
in x- and y-directions, things can be different, as will be 
discussed in the next section. 

External magnetic field will shift the chemical poten- 
tial of the component m into fi m = /j,q + mEz, where Ez 
is the Zeeman shift and is the chemical potential of the 
?7i = component. An external magnetic field will tend 
to magnetize the system and move more population into 
the state with the lowest chemical potential. It is simple 
to show, that for a non-rotating homogeneous spinor con- 
densate the Zeeman shift required for a complete magne- 
tization in the direction of the external magnetic field is 
on the order of \ a n, where n is the density p3]. Note that 
the magnetizing field is very small e.g. for sodium atoms 
with reasonable densities it is about B = 10 -4 — 1CP 3 G. 

Assume a slowly rotating system prepared in the ab- 
sence of the magnetic field and allowed to equilibriate 
into the ground state. We then turn on the magnetic 
field which is aligned with respect to the rotation axis 
and the angular momentum of the system. With an in- 
creasing magnetic field the ground state of the system 
evolves towards the ground state of the scalar conden- 
sate in m = 1 or m = — 1 component. Which component 
is favored depends on the direction of the magnetic field. 
As the rotation frequency is small, this ground state has 
a vanishing angular momentum. Therefore angular mo- 
mentum of the ground state decreases with increasing 
magnetic field. 

The reverse process of starting with a large magnetic 
field along the z-direction and then turning it off behaves 
differently. The ground state of the scalar GP equation 
is metastable under such change. Initially with a large 
magnetic field the ground state is essentially a scalar con- 
densate without angular momentum. It is hard to relax 
such a state into a coreless vortex state with finite an- 
gular momentum. This behavior is reminescent of the 
hysteresis for the vortex nucleation in scalar condensates. 

VI. CREATION 

As it was already mentioned, the behaviour of the sys- 
tem may be different if the external magnetic field has 



other non-zero components than the one pointing in z- 
direction. This is dramatically shown in the topologi- 
cal vortex formation [B4 Hj7| which was recently demon- 
strated experimentally P8|| T In this experiment a mag- 
netic field configuration of the Ioffc-Pritchard trap was 
changed in time. In particular the z-component of the 
magnetic field was inverted. Once the z-component is in- 
verted, the condensate wavefunction has acquired a wind- 
ing number 2. It should be noted that this process con- 
tinuously transforms a condensate with vanishing angular 
momenta into a condensate with an angular momentum 
2 per particle. When the spreading of the wavefunctions 
can be ignored this process is equivalent with coreless 
vortex transformations. In this respect coreless vortices 
have already been experimentally created. In a Ioffc- 
Pritchard magnetic trap a coreless vortex is not stable 
since all hyperfine states are not trapped. The stabiliza- 
tion of the coreless vortex structures would require an 
optical trap. 

VII. CONCLUSIONS 

We have studied the ground state properties of rotat- 
ing spin-1 condensates. We have found out that for a 
sufficiently small rotating frequency the ground state of 
a ferromagnetic spinor condensate is a coreless vortex. 
However, because of the finite size of a condensate a cer- 
tain critical rotation frequency has to be exceeded for 
a coreless vortex to appear. We have also investigated 
antiferromagnetic condensates which in contrast to core- 
less vortices have rotating ground states exhibiting the 
discrete staircase behavior. 

All the simulations presented in this paper have been 
performed in two dimensions. Consequently they are 
expected to describe well pancake shaped condensates. 
We have also performed simulations with full three- 
dimensional Gross-Pitaevskii equations and convinced 
ourselves that our results are not specific to two dimen- 
sions. While a detailed study in three dimensions is at 
present too time-consuming, it does seem that our results 
apply qualitatively even in a cigar shaped trap. 

In this paper we have focused solely on the ground 
state. Naturally the ground state is one thing, but get- 
ting the system into it in a reasonable time is quite an- 
other. The energy landscape of the rotating spinor con- 
densate is complex and metastable states abound. This 
may cause problems in the experiments as the relaxation 
time to the true ground state can become very long. Also 
the problem of nucleation of coreless vortices may cause 
problems, especially if the system is initially prepared as 
a pure stationary condensate. Nevertheless with some 
effort all these complications can be overcome, and the 
ground state of the rotating spinor condensate promises 
to be an exciting field of research. 
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